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Improved QCD expansion and the problem of determination of and 
the condensates from the r decays 
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The problem of determination of a, and the condensates from the moments of invariant mass distribution in the 
semileptonic r decays is considered. It is investigated how the extracted values of Qf3(m^) and the condensates 
are affected by the the renormalization scheme dependence of the next-to-next-to-leading order perturbative 
corrections to the spectral moments. A simplified approach is used, in which the nonperturbative contributions 
are approximated by the terms of dimension six in the SVZ expansion, which arise from the four-quark condensates. 



Recently there has been considerable interest in 
the determination of the QCD parameters from 
the invariant mass distribution in the semilep- 
tonic r decays . It appears that using specific 
moments of the invariant mass distribution one 
may obtain rather tight experimental constraints 
for Us and the QCD condensates. The perturba- 
tive QCD corrections to the spectral moments — 
which are known to next-to-next-to-leading order 
(NNLO) — play central role in such a determi- 
nation. These corrections are usually evaluated 
in the MS renormalization scheme. However, in 
NNLO there is a two-parameter freedom in the 
choice of the renormalization scheme (RS). The 
numerical values of the perturbative corrections 
depend on the choice of the scheme. It is there- 
fore an interesting question, how the values of 
the QCD parameters, extracted from the r decay, 
may be affected, if a different choice of the RS is 
made, motivated for example by the principle of 
minimal sensitivity (PMS) ||^. 

The theoretical framework adopted in the anal- 
ysis of the T decay data involves the -R^V/A 
ments, defined by the relation 
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where Fg is the electronic width of the r lep- 
ton, dT^li^ /ds is the invariant mass distribution 
for the Cabbibo allowed semileptonic t decays in 
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the vector (V) or axial-vector (A) channel and 
w^^{s) = (1 - s/ml)''{s/miy. The QCD predic- 
tions for Rfy/ji^ have the form: 
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where = 0.9752 and Sew = 1.0194. The 

i?Q'' factor denotes the parton model prediction. 

The dpi term is the perturbative contribution, 

which is evaluated assuming 3 massless quarks — 

in this approximation it is universal for the V and 

A channels. The (5*^' term denotes the Con- 
ner y I A 

tribution from the nonperturbative QCD effects, 
which are estimated using the SVZ approach 
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D=4.6.. 



where < O^^^ > are the vacuum expectation val- 
ues of the gauge invariant operators of dimension 
D and c^' are coefficients specific for the consid- 
ered spectral moment and the type of the opera- 
tor. 

The authors of |,| used ^nd the -R^J^/^ 

moments with k — 1 and 1 — 0,1, 2, 3, and fitted 
the Z? = 4, 6, 8 condensates. Since our goal is pri- 
marily to study the theoretical uncertainties in 
the whole procedure, we shall adopt a simplified 
approach, in which the nonperturbative contri- 
bution to R^\/j^ is approximated by the lead- 
ing D — 6 contribution, coming from the four- 
quark condensates. (This means we neglect the 
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D = A contribution, which is suppressed by an 
additional power of as, the higher order pertur- 
bative corrections to the coefficients of the D = 6 
operators, and the D > 8 contributions.) We 
also use the R].\/j^ moment, for which a simi- 
lar approximation may be made, and for which 
we have the relation <5(|)y/^ = —(70/13)5™^^^^. 
This gives us two equations to fit two parame- 
ters: A^IL (or, more conventionally, af^'^(m^)) 
and one parameter characterizing the nonpertur- 
bative contribution, which we take simply to be 

°(6)yM- 

The perturbative QCD corrections S^l are eval- 
uated using a contour integral expression, which 
relates them to the QCD correction to the so 
called Adler function, i.e. the logarithmic deriva- 
tive of the transverse part of the vector/axial- 
vector current correlator: 



^.1 - - I -f' (4) ^^*(-). (4) 

where f''\a/m'^) is a weight function specific 
to the considered moment and C in our case 
is assumed to be a circle a = — exp(— i6'), 
d € [— TT, tt]. The NNLO renormalization group 
improved perturbative expansion for 6^ may be 
written in the form: 

= a{-G)\\ + ria(-a) -|- r^a'{~a\ (5) 

where a = asji: = /(Ai:'^) denotes the running 
coupling constant that satisfies the NNLO renor- 
malization group equation: 

da b 

By evaluating numerically the integral (Q) with 
the renormalization group improved expression 
for 5^ under the integral we resum to all orders 
some of the corrections, which would appear in 
the "naive" expansion of in powers of as{rn^) 

i- 

The coefficients ri, r2 and C2 are RS dependent, 
but there exists a RS invariant combination: 



{1 + Cia + C2a'^). 
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/32 = C2 + r2 - ciri 



(7) 



For the 5n we have fSp = 5.23783. We shall 
use ri and C2 to parametrize the freedom of choice 
of the RS in NNLO. 



The dependence of 5^ on the scheme parame- 
ters ri and C2 was discussed in detail in [|lo| and 
the RS dependence of 5^ was investigated in . 
In both cases it was found that for moderate val- 
ues of ^jjg the NNLO predictions have a saddle 
point type of behavior as a function of ri and C2 
and that the position of the saddle point is well 
approximated by ri = and C2 = 1.5y02 = 7.85 7. 
For very large values of Aii^. the RS-dependence 
pattern is more complicated than a simple sad- 
dle point, but even then the scheme parameters 
distinguished above belong to the region of ex- 
tremely small RS dependence. We shall therefore 
accept these parameters as the PMS parameters 
in NNLO. 

In order to obtain a quantitative estimate of 
the RS dependence of the results we shall use the 
approach outlined in | [T^ , based on the existence 
of the RS invariant p2- In ||l2| it was proposed 
to calculate variation of the predictions over the 
set of schemes for which the expansion coefficients 
satisfy the condition: 

(T2(?-l,r2,C2) < I \p2\, (8) 

where 

CT2(ri,r2,C2) = |c2| -t- \r2\ +ci|ri| + rl. (9) 

A motivation for the condition (||) is that it elimi- 
nates schemes in which the expansions 

(|) and (§) 

involve unnaturally large expansion coefficients, 
that introduce large cancellations in the expres- 
sion for the RS invariant p2- The constant I in 
the condition (||) controls the degree of cancella- 
tion that we want to allow in the expression for 
P2. In our case we have for the PMS parameters 
(T2(PMS) « 2|p2|, so in order to take into account 
the schemes, which have the same — or smaller 
— degree of cancellation as the PMS scheme we 
take I — 2. 

Let us begin with the fits in the vector chan- 
nel. We use the experimental values reported 
recently by ALEPH g: R^^y = 1.782 ± 0.018, 
Dl\ = Rl^y/R"J^y = 0.0532 ± 0.0007. For sim- 
plicity we neglect correlations between experi- 
mental errors for and D^^. For comparison 
with other determinations of ag we evolve the fit- 
ted value of a^^{m'^) to the scale of using the 



NNLO renormalization group equation and the 
NNLO matching formula [|l3| at = 2mc, 2mi,. 

Using the PMS predictions we obtain from the 
NNLO fit in the vector channel S"^°^y = O.Of 56 ± 

0.0023 and = 42f ± 30MeV, which corre- 

MS ' 

sponds to as(TO^) = 0.356 ± 0.017 and as{ni^) — 
0.1226 ± 0.0018. For comparison, using the MS 
scheme we obtain = 441 ± 32 MeV and 

MS 

joOj^ = 0.0147 ± 0.0025, which corresponds to 

as{ml) = 0.367 ± 0.018 and as{ml) = 0.1238 ± 
0.0018. 

In order to make our calculations more gener- 
ally useful we show in Fig. j| the results of the 
NNLO fit of as(TO^) and S'^-^y, obtained using 
the PMS predictions, as a function of the experi- 
mental values of i?^^ and D].\. 
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In Fig. H we show how the value of asim^) 
resulting from the fit depends on the param- 
eters ri and C2 specifying the renormalization 
scheme in NNLO. By varying the RS param- 
eters in the region satisfying the condition (|^) 
with ^ = 2 we obtain the variation 0.347 < 
asiml) < 0.367 (403 MeV < < 440 MeV, 

0.1217 < as{ml) < 0.1238). Analyzing in a sim- 
ilar way the RS dependence of S'^^^y we obtain 
0.0145 < S"° < 0.0182. 
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Figure 2. The contour plot of the fitted value 
of as(w^) in the vector channel as a function 
of the RS parameters ri and C2- The region of 
the scheme parameters satisfying the condition 
) with I — 2 has been also indicated. 



Figure 1. Plot of the fitted values of Q!s(m^) and 
(5°g^ in the vector channel as a function of R'^^y 

and Dl'^y, obtained using the NNLO PMS pre- 
dictions. The dashed lines indicate the change in 
the plot when the MS NNLO predictions are used 
instead. 



It is of some interest to perform the same fits 
using instead the NLO predictions. (The PMS 
parameters in NLO are ri = —0.76 for (5°^ and 
ri = —0.64 for 6^^.) Using the NLO predictions 
in the PMS scheme we obtain ,5™ ^ = 0.0150 and 

A^ = 465 MeV. Using the NLO predictions in 
the MS scheme, we obtain (59^^^ = 0.0148 and 



AiiL = 527MeV. 

A'/ S 

Similar fits may be performed in the axial- 
vector channel. We use the experimental val- 
ues reported recently by ALEPH @: = 
1.711 ± 0.019, 1)12^ = 0.0639 ± 0.0005. Using 
the PMS scheme we obtain in NNLO S^^^-^^ = 

'^^^ - 380 ± 34MeV, 
0.335 ± 0.018 
0.1203 ± 0.0021. For compar- 



-0.0165 ± 0.0018 and 
which corresponds to as(rn^ 



and as{rn^) 

ison, using the MS scheme we obtain (5?g, 



-0.0168 ± 0.0021 and = 398 ± 37MeV, 

M S 

which corresponds to asim^) = 0.344±0.019 and 
a,(m|) = 0.1213 ±0.0021. 

Performing the variation of the RS parame- 
ters in the region satisfying the condition (^) 
with ^ = 2 we obtain for the axial-vector chan- 
nel 0.326 < as{ml) < 0.343 (364MeV < aI^L < 
397 MeV, 0.1193 < as{ml) < 0.1212) and 
0.015 < -(5['gO ^ < 0.017. 

Summarizing, we find that in NNLO the change 
from the MS scheme to the PMS scheme results in 
the reduction of the fitted value of as(m^) by ap- 
proximately 0.01. (as(m|) is reduced by 0.001.) 
Also, the difference between the NLO and NNLO 
results is much smaller in the PMS scheme than 
in the MS scheme. The extracted values of the 
QCD parameters appear to be relatively stable 
with respect to change of the RS. However, vary- 
ing the scheme parameters ri and C2 in the region 
satisfying the condition (||) with Z = 2 we obtain 
an uncertainty in the extracted value of as{'m'^) 
of approximately 0.02 (uncertainty in as(m\) is 
0.002). 

More details about the calculations and the re- 
sults reported above may be found in p[. 
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